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Abstract 

Background: Niche construction has received increasing attention in recent years as a vital force in evolution and 
examples of niche construction have been identified in a wide variety of taxa, but viruses are conspicuously absent. 
In this study we explore how niche construction can lead to viruses engineering their hosts (including behavioural 
manipulation) with feedback on selective pressures for viral transmission and virulence. To illustrate this concept we 
focus on Baculoviridae, a family of invertebrate viruses that have evolved to modify the feeding behaviour of their 
lepidopteran hosts and liquefy their cadavers as part of the course of infection. 

Results: We present a mathematical model showing how niche construction leads to feedback from the behavioural 
manipulation to the liquefaction of the host, linking the evolution of both of these traits, and show how this 
association arises from the action of niche construction. Model results show that niche construction is plausible in this 
system and delineates the conditions under which niche construction will occur. Niche construction in this system is 
also shown to be sensitive to parameter values that reflect ecological forces. 

Conclusions: Our model demonstrates that niche construction can be a potent force in viral evolution and can lead 
to the acquisition and maintenance of the behavioural manipulation and liquefaction traits in Baculoviridae via the 
niche constructing effects on the host. These results show the potential for niche construction theory to provide new 
insights into viral evolution. 



Background 

Niche construction, defined as the "process whereby 
organisms, through their metabolism, their activities, 
and their choices, modify their own and/or each others 
niches" [1], has received growing attention as a potentially 
potent force in evolution. This act of modifying niches 
changes biotic and abiotic sources of selection in the 
organismss environment, introducing feedback into the 
evolutionary process as the effects of niche construction 
modify an organisms genetic and ecological inheritance 
[2,3]. Examples of niche construction have been identified 
in all taxonomic groups, but little attention has been paid 
to viruses. As obligate intracellular parasites, viruses nec- 
essarily have ongoing (and often antagonistic) interactions 
with their hosts; by altering their hosts they are altering 
their own selective environments. When these interac- 
tions modify the host in a way that has lasting effects 
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on viral selection pressure, niche construction will be in 
operation. 

In a fashion similar to other parasites [4], host manip- 
ulations that create niche construction effects for viruses 
may occur at multiple levels of the host organisms. For 
example, at the level of the hosts physiology, viruses 
such as norovirus cause fever, diarrhoea, and vomiting. 
These can drive niche construction because of their direct 
effect on the host (e.g. transmission via viral shedding in 
diarrhoea), and there may be additional effects from the 
hosts response to these changes. Adelman and Martin [5] 
describe the behavioural changes in the host from fever 
and other symptoms as an emergency life-history stage in 
which organisms prioritise recovery from infection; the 
host response, such as lethargy or anorexia, may mod- 
ify the virus' own selective pressures. Viruses also directly 
cause gross changes in host behaviour, such as induction 
of pathological aggression in rabies, or reports of pro- 
longed feeding by mosquitoes infected with dengue virus 
[6]. Viruses even have lasting effects on the environment 
beyond the host, as when marine viruses affect oceanic 
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participation in carbon cycling [7]. Of course, this process 
of modifying the host may also change the hosts selective 
pressures; beyond the response of the immune system, the 
host may be under selection to counter the effects of the 
physiological or behavioural manipulations of the virus. 
Our understanding of viral evolution will be enhanced 
by identifying and examining traits that are coupled in 
this way, especially traits through which niche construc- 
tion generates feedback into viral life history properties. 
In addition, the study of niche construction as an evo- 
lutionary process will be enhanced by new avenues for 
empirical tests in viruses. For example, we might com- 
pare populations of viruses in which the hypothesised 
niche construction pathway (such as a host manipulation 
enhancing transmission) is operating with populations in 
which this pathway has been experimentally negated to 
compare fitness or the evolution of viral traits [3]. With 
their short generation time and quick rates of evolution, 
viruses are attractive targets for empirical tests of niche 
construction theory. 

In this study, we illustrate the operation of viral niche 
construction with an example drawn from the family 
of invertebrate viruses known as Baculoviridae. Bac- 
uloviruses are rod-shaped, enveloped dsDNA viruses with 
large genomes that range in size from 90 kb to 180 
kb [8]. Two major baculovirus morphologies are dis- 
tinguished, the nucleopolyhedroviruses (NPVs) and the 
granuloviruses (G Vs) based on the contents of their occlu- 
sion bodies (OBs, the proteinaceous structures that con- 
tain the virions); NPVs have multiple virions per OB, 
while GVs typically have only one. NPVs are further sub- 
divided into single or multiple NPVs (SNPV / MNPV) 
based on the contents of their virions, which can hold 
a single or multiple nucleocapsid. Baculoviruses infect 
arthropod hosts and have been reported worldwide from 
over 600 host species, primarily from the order Lepi- 
doptera though they have been also identified in Diptera, 
Hymenoptera, and Decapoda as well [8]. This host speci- 
ficity forms the basis for the currently accepted classifi- 
cation of the family Baculoviridae into four genera [9]: 
Alphabaculovirus, lepidopteran-specific NPVs, Betabac- 
ulovirus, lepidopteran-specific GVs, Gammabaculovirus, 
the hymenopteran-specific NPVs, and Deltabaculovirus, 
the dipteran-specific NPVs . Baculoviruses are evolution- 
arily old, with molecular estimates placing their origin 
around 310 million years ago during the Carboniferous 
period of the Paleozoic era [10]. 

In many lepidopteran host species it has been known 
since at least the 19th century that baculoviruses 
induce a behavioural modification in their hosts which 
results in forced climbing behaviour (referred to as 
"Wipfelkrankheit" or tree-top disease; Hoffman, cited by 
Goulson [11]). As part of their pathology, baculoviruses 
are also known to liquefy their hosts after death, allowing 



occluded virions to burst and enhance their dispersal [12]. 
Mechanistically, a small set of viral genes are known to 
be directly involved in these actions. Recent work has 
shown [13] that compared to controls, deleting the viral 
gene ecdysteroid UDP-glocosyltransferase(egt) in LdM- 
NPV (host: Ly mantra dispar) will disrupt the climbing 
behaviour and reinstating the gene also restores the climb- 
ing behaviour. Other genes, such as the viral gene for 
the baculovirus protein tyrosine phosphatase (ptp) of 
BmNPV (host: Bombyx mori), are occasionally implicated 
in the behavioural manipulation as well [14]. Important 
viral genes for liquefaction are chitinase and cathepsin 
(e.g. [15,16]), though other candidates may also be 
involved [17]. Genes such as egt, chitinase, and cathep- 
sin are not necessarily the only causes for these traits, 
and given the constraints on viral genomes it is unsurpris- 
ing that they may also have other functions (e.g. egt and 
virulence: [18,19]). 

Despite the advances in understanding the mech- 
anistic basis for host manipulation and liquefaction, 
insight into the evolutionary consequences and origins 
of these baculovirus characteristics remains speculative. 
It is thought that tree-top disease combined with liq- 
uefaction enhances transmission of the virus to con- 
specifics by shedding virus directly into the path of feeding 
conspecifics (as the virus infects orally). However, this 
leaves open the question of how such a conjunction of 
behaviours could arise to begin with. As Hoover et al. [13] 
note about L. dispar, in the absence of the virus cater- 
pillars will hide in bark crevices or climb down the tree 
during the day to avoid predation from birds; this sug- 
gests that climbing behaviours without liquefaction would 
harm the transmission of the virus. Similarly, liquefaction 
of the host may have a much greater effect on transmission 
if it occurs at the top of the host plant / tree, as the virus is 
spread into the path of feeding conspecific by mechanical 
factors such as wind and rain. In contrast, liquefaction at 
the bottom of the host vegetation or in crevices may lead 
to sequestration of the virus in the soil and bark reservoir 
that leads to next seasons infections but reduce trans- 
mission in the current season, as virus left unprotected 
will degrade quickly from UV exposure or be removed by 
other mechanical means [20]. Thus, each feature of bac- 
ulovirus pathology may be slightly beneficial, neutral, or 
even deleterious on its own, but feedback from one feature 
to the other may link the evolution of the two together. 

In this study we aim to contribute to the broadening of 
the domain of niche construction to viral evolution. To do 
so, we apply a niche construction perspective to make a 
general evolutionary argument about the acquisition and 
maintenance of tree-top disease and liquefaction in bac- 
uloviruses infecting Lepidopteran hosts. We construct an 
epizootic model of baculovirus transmission to link viral 
fitness to genotypes with each combination of features and 
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demonstrate that there are valid evolutionary paths from 
a hypothetical virus with neither feature to one such as 
LdMNPV with both (Figure 1). 

Methods 

Mathematical model of viral niche construction 
Epizootic model 

To study the effects of baculovirus behavioural manipu- 
lation and liquefaction on viral fitness, we begin with a 
compartmental epizootic model similar to a Susceptible- 
Exposed-Infectious-Recovered (SEIR) model [21]. Our 
model tracks the status of the host population as it transi- 
tions through five states (c.f. [20,22]): susceptible individ- 
uals (S) that have not been exposed to the virus, exposed 
individuals (E) that have ingested the virus, intact cadav- 
ers (C) that spread the virus at a low rate, liquefied cadav- 
ers (L) that spread virus at a higher rate, and a reservoir 
(V) that harbours virus for the next epizootic outbreak 
(tracked as proportions, not absolute numbers of individ- 
uals). Epizootics are assumed to begin when some small 



proportion of the recently hatched larvae come in con- 
tact with virus from the previous seasons reservoir and 
become infected (E(0) ^ 1.0 x 10 -5 ). Susceptible indi- 
viduals enter the exposed class by coming in contact with 
intact or liquefied cadavers, after which they die at some 
rate v. 

Intact cadavers infect susceptibles at a rate /3c, are liq- 
uefied at a rate k, and shed virus directly into the reservoir 
at a rate 0. The egt gene also increases the virion yield by 
deactivating moulting hormone and causing infected lar- 
vae to grow larger and die later than uninfected larvae. 
This potentially niche constructing effect can be incorpo- 
rated into the transmission parameters /3l and /3c, as an 
increased virion yield would increase transmission. 

As we are tracking the infection for a single season, we 
presume a host in which eggs overwinter and hatch in 
spring; thus, we neglect births in the model, and track 
only baseline mortality from natural causes and preda- 
tion (/x). However, we do allow for excess mortality as 
a consequence of viral infection. For instance, exposed 
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Figure 1 Schematic representation of niche construction. Conceptual representation of the action niche of construction in the baculovirus 
system. The build up of EGT within a lepidopteran host over successive viral generations (A) is an instance of ecological inheritance, which leads to a 
change in the host behaviour by time T (the 'zombie' trait of the virus; see Methods). This enhances the transmission of the virus if it liquefies the host 
(the 'gooey' trait). As described in 'Genetic structure' below, we study four viral genotypes encoding these two traits: zg (non-zombie, non-gooey), 
Zg (zombie, non-gooey), zG (etc.), ZG; the genotype fitnesses in this figure are arbitrary chosen for illustrative purposes. In (B) we represent a 
potential pathway through genotype space for the virus over evolutionary time. (A) shows the change in fitness landscape as a result of the viral 
niche construction (i.e. its host manipulation), while (B) shows the evolutionary consequences of this niche construction activity for the virus. 
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hosts suffer from reduced mobility and behavioural 
changes that may increase their chance of being predated 
(e.g. [13,23]) or cannibalised [24], and so this rate which 
we denote 8 is higher as the disease progresses (8 > /jl). 
Intact cadavers that die are liquefied internally until their 
cuticle bursts in a matter of hours; during this time, they 
can be removed by predation [12] at the same 8 rate as 
intact cadavers, after which they transmit to susceptibles 
at a rate fii, and contribute to the reservoir at a rate r. 
Virus that enters the reservoir degrades at some small rate 
6. Note that this model includes a protected reservoir that 
will cause future outbreaks sensu [22] but, as shown by 
[25], simple environmental persistence of the virus could 
be used in the model instead with the same result. A flow 
diagram of this model is given in Figure 2. We focus on 
the dynamics of acute infection by baculovirus of a local 
population of larval hosts and consider long-term niche 
construction arising from these dynamics. 

We model the movement of hosts between compart- 
ments as a set of ordinary differential equations as follows. 



S = -foLS - p c CS 

E = p L LS + p c CS - Qi + v)E 

C = vE-(k + 6 + 8)C 

L = kC - (t + 8)L 

V = 0C+rL-eV. 



(1) 
(2) 
(3) 
(4) 
(5) 



The parameters for this model and their meanings can 
be found in Table 1 (parameters values are discussed 
below). We follow [21] and [26] (see also [27]) to calculate 
the basic reproduction number Ro of this model using the 
next generation matrix, as follows. 

We order the compartments of the model as 
(£, C, L, S, V) so that the first m = 3 contain infected indi- 
viduals. Let x = (x±, . . . ,x n y,Xi > 0 be the corresponding 



number of individuals in each disease compartment. Let 
J~i(x) be the rate of new infections appearing in compart- 
ment i, and Vi(x) be the net transfer of individuals into 
compartment /. The difference Ti(x) — Vi(x) gives the 
rate of change of X{. 
In our model, 



(p L LS + p c CS\ 
0 
0 
0 
0 



and{V/} = 



/ Qi + v)E \ 
■vE+(k+0+8)C 
-kC+(r+8)L 

\-6C -zL + eV ) 



(6) 



We form m x m matrices F and V whose elements are 
given by 



hi 



~dFi(x 0 y 


and Vy = 




dxj 




dXj 



(7) 



where xq is the disease-free equilibrium (here, xq = 
[0,0,0,1,0]), and/,; < m. 
Using equation (7), F and V are 

/0 Pc Pl\ 
F= 0 0 0 
\0 0 0 / 



and 



/fJL + V 0 0 \ 

V = I -v k + 0+8 0 

\ 0 —X r + aJ 



The basic reproduction number Rq of the model is the 
spectral radius of the next generation matrix which is 
given by FV~ l , which simplifies to 

p c v(8 + x) + p L v\ 



Ro 



(5 + ^+A)(/x + v)(5 + r). 



(8) 



V 5 5 




e 



Figure 2 Model diagram. Model diagram. Susceptibles are infected at rates /3c and fii from intact and liquefied cadavers respectively. Exposed 
hosts die from infection at rate v and from natural causes or predation at rate /x. Intact cadavers are liquefied at rate k and contribute to the 
reservoir at rate 9. Both intact and liquefied cadavers are removed by predators and other causes at the same rate 8. Liquefied cadavers contribute 
to the reservoir at rate 6, and virus in the reservoir degrades at rate e. The parameters are summarised in Table 1 . 
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Table 1 Parameter values 



Parameter Meaning 

Pi Transmission between susceptibles and liquefied cadavers. 

Pc Transmission between susceptibles and intact cadavers. 

Ii Baseline death rate for susceptible and exposed 

(but asymptomatic) hosts. 
v Rate of conversion from exposed to intact cadaver states. 

8 Clearance rate for intact and liquefied cadavers. 

Includes predation and cannibalism that does not transmit 

infection. 

X Rate of conversion of intact to liquefied cadavers, 

r Rate at which virus shed from liquefied cadavers enters 

the reservoir. 

6 Rate of virus shedding into the reservoir from intact 

cadavers. 

6 Clearance rate of virions from the reservoir by 

environmental factors (UV degradation, wind, rain, etc). 



We can clarify the contribution of each model component 
by rewriting this expression as 



Pc 



+ 



8 + e+X\v + nJ S + t \v + nj \8 + 9+X 



(9) 



The first term above gives the number of new infec- 
tions produced by an intact cadaver, over (8 + 0 + A.) -1 
time units. It includes the factor v/(v + fi) which is the 
proportion of exposed individuals that die. The second 
component corresponds to the contribution of the lique- 
fied state and includes both the proportion entering the 
cadaver state from the exposed state and the proportion 
entering the liquefied state from the cadaver state, given 
byA/(A + <5+0). 

Genetic structure 

To trace the selective effects of the behavioural manipu- 
lation and liquefaction features of baculovirus, we define 
four genotypes from two hypothetical loci with two alleles 
each. For the purposes of our model we adopt a simplified 
view of egt and chitinase / cathepsin as being responsible 
for the host manipulation and liquefaction traits respec- 
tively (Figure 1); further, we adopt non-specific names for 
these traits to highlight the general nature of our evo- 
lutionary argument which shows the niche constructive 
pathway from a historical but unknown protovirus with 
neither trait to a modern virus with both. Z and z are zom- 
bie and non-zombie alleles, corresponding to the presence 
or absence of the behavioural manipulation driven by the 
egt (and ptp) gene. G and g are gooey and non-gooey alle- 
les, corresponding to the presence of the liquefying genes 



cathepsin and chitinase. This gives us zg Zg zG, ZG as 
the set of genotypes in the system, each with its own basic 
reproduction number Roj where i denotes the genotype 
(with i e {zg, Zg, zG, ZG}), 

As shown in Table 2, we define the effect of the geno- 
types by their effect on viral life history traits via model 
parameters, with each genotype having its own value of 
four critical parameters indexed by genotype: namely, for 
the i th genotype, fin, T;, A.;, and 0,-. Note that we replace the 
generic symbols introduced in the section above describ- 
ing the ecology of disease spread with new notation given 
in Table 2, which depend on genotype. Here we introduce 
a new parameter, y, which represents the advantage or 
disadvantage to each mode of transmission (horizontal or 
via the reservoir in next season). Relative to the zg geno- 
type ZG increases horizontal transmission by liquefying 
'tree-top' (plant-top, etc) larvae and increasing spread of 
occluded virus over foliage that conspecific larvae will eat 
(increasing Pl,zg)- However, this reduces contribution to 
the reservoir (lowering tzg) because exposed virus will 
degrade quickly, especially as a result of UV light ([28], 
p.364-5); in the past, this fact limited the potential use of 
baculoviruses as insecticides. We assume that the reverse 
happens to the zG genotype, in which insects die in a nat- 
ural position closer to the soil or understory and then 
liquefy to spread virus to the reservoir (soil, understory, 
tree bark, etc). Occluded NPVs in forest settings have 
been shown to persist as long as 41 years after a natu- 
ral outbreak [29]. Thus, y decreases Pl,zG and reduces 
t z g relative to the zg genotype. Since they do not liquefy 
their host (k zg , kzg = 0), it is assumed that the Zg and zg 
genotypes are neutral with respect to transmission of the 
virus (relative to the protovirus zg). However, the zombie 
allele would still modify contribution to the reservoir by 
the same argument as above, so we can represent this by 
modifying the values of 0/; we represent this as a low value 
of 0/, 6 1 for the Z allele and a high value 0^ for the z allele 
{0 h > 0/). 

Based on the principle of competitive exclusion in 
pathogen virulence [30], we assume that evolution fol- 
lows a path in which viruses of higher Rqj invade and 
replace those of lower 7?o,/; this approach is widely used 
in evolutionary epidemiology, as noted by [31]. Com- 
petitive exclusion also allows us to ignore the dynamics 



Table 2 Genotypes 


/ fiu 






Oi 


ZG Pl0+y) 


r(1-y) 


X 


Oi 


Zg p L 


T 


0 


01 


zG p L Q -y) 


r(1+K) 


X 




zg p L 


T 


0 





Genotypes and their effect on the model parameters. 
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of the reservoir, as we assume that the dominant geno- 
type will also be dominant in the reservoir, and the sea- 
sonal dynamics of baculovirus have already been modelled 
(e.g. [22]). Thus, when the ZG genotype has the highest fit- 
ness value, we conclude that the trajectory through geno- 
type space will proceed according to the changed selective 
pressures imposed by niche construction (Figure 1). If 
ZG does not have the highest fitness in some region of 
the parameter space, we conclude that niche construction 
does not occur. Thus, using values of the basic reproduc- 
tion number 7?o,/ as a function of genotype, we explored 
the potential for niche construction in parameter space by 
defining fitness orderings by genotype. Plotting viral fit- 
ness over y allows us to rank the value of Roj for each 
genotype at a given set of parameter values in order from 
lowest to highest, yielding a fitness ordering to which we 
assign a value and a colour (Table 3). 

Parameter values 

For numerical analyses of the model, we set the baseline 
parameters as described in Table 1 as follows (all units are 
rates per individual per day). Parameters /3l and /3 c were 
set to 1 and 0.5 per day respectively, to provide a large 
transmission advantage to liquefied cadavers, following 
baculovirus pathology and ecology [12]. Contribution to 
the reservoir, r (set to 0.2) and 0 (set to 0.25), is expected 
to be low (conditional on where the virus is shed). For 
the remainder of this study we neglect mortality to sus- 
ceptible and exposed hosts (//,) and clearance from the 
reservoir (e) as being unimportant for the niche con- 
struction aspect of the model and so set them to zero, 
but deaths from predation may have an effect on viral 
fitness and so we set a low value of 8 = 0.05 as a baseline. 
Finally, we take the average time from exposure to death 
as being 5-10 days (as is typical for fourth- or fifth-instar 
infections [23]), and so set the rate of disease progression 
from exposed to intact cadaver (v) to 1/7 = 0.14. In 
exploring the effect of model parameters on niche con- 
struction, we varied fii from 0.5 to 2, /3c from 0.1 to 1, 0/ 
from 0.01 to 0.5, X from 0.1 to 2, 8 from 0 to 0.5, and r 
from 0.01 to 0.5. 



Table 3 Color mapping 

Fitness ordering Color 

Any in which ZG doesn't have the highest fitness. Purple 

zG < Zg < zg < ZG Blue 

Zg < zG < zg < ZG Teal 

Zg < zg < zG < ZG Green 

zG < zg < Zg < ZG Yellow 

zg < Zg < zG < ZG Orange 

zg < zG < Zg < ZG Red 



Mapping of the niche construction evolutionary pathways to colors. 



Model analysis 

We numerically solved equations (l)-(5) using Mathe- 
matica [32]. Initial conditions were set to C(0) —> 1 x 
10- 8 ,£(0) -> 0,1(0) -> 0,V(0) -> 0, andS(0) -> 1- 
(£(0) + C(0) + 1(0) + V(0)). Calculation and plotting 
of genotype fitness and fitness orderings over parameter 
space was done using R [33]. 

Results 

Model results 

Figure 3 shows the course of infection under the model 
using the parameters values given in Section parame- 
ter values! These dynamics are qualitatively in line with 
known properties of baculovirus infections of species such 
as L. dispar (e.g. [34,35]). Larvae of L. dispar hatch in early 
spring and proceed through five or six instars of roughly 
two weeks each [36]. 

Viral fitness (Rq) as a function of genotype and the 
advantage parameter y is shown in Figure 4. For this com- 
bination of parameters, low values of y (less than 0.3) lead 
to evolution proceeding from the protovirus with neither 
trait (zg) and acquiring either trait first before arriving at 
the zombie-gooey genotype ZG via the niche constructive 
feedback on the virus's selective pressures. At higher val- 
ues of y (greater than 0.5), we see a fitness ordering that 
only allows niche construction to proceed by acquiring the 
zombie trait first (ZG). 

Figure 5 shows the plausibility of niche construction 
under various model parameters. Our results show that 
niche construction is plausible in this system for a wide 
range of parameters; we also see that some combina- 
tions of parameters do not allow niche construction. The 
rates of transmission (fii, f3c), combined with the rate 
at which the virus is shed into the reservoir (r, 0) from 
either source affect the sequence of trait acquisition that 
is most profitable for the virus, such that low / high trans- 
mission from liquefied and intact cadavers respectively 
restricts the virus to acquiring the zombie trait to increase 
its fitness when the contribution to the reservoir is high 
(Figure 5(a,c)). The value of 0/ shows the relative contribu- 
tion to the reservoir from the Z allele compared to the z 
allele. Because 0^ was set to 0.25 (Table 1), the vertical axis 
of Figure 5(b) is the difference between the Z and z alleles. 
Thus, when 0[ goes above 0.25, the effect of the zom- 
bie trait inverts. The interplay between these parameters 
defines regions in which the transmission from lique- 
fied cadavers shows a trade-off with the contribution to 
the reservoir made by zombie caterpillars. When 0/ is 
low (large effect of Z) and transmission is low, no niche 
construction occurs, but high transmission rescues the 
niche constructive pathway. When the effect of Z inverts, 
high rates of transmission still allow niche construction 
to acquire both traits (green region), but only through 
the gooey trait. And at low rates of transmission, niche 
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Figure 3 Epizootic dynamics. Epizootic dynamics over time, showing the progress of the infection in a susceptible population that has been 


exposed to a single infectious individual. For this figure, ft = 1 , ft = 0.5, z = 0.2, 6 = 0.25, v = 0.14, 8 = 0.05, [i = 0, e = 0. 
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Figure 4 Viral fitness as a function of y. Plot of viral fitness (measured by R 0 ) as a function of the advantage parameter y. For this figure, ft = 1 , 


ft = 0.5, r = 0.2, 6 = 0.25, v = 0.14, 8 = 0.05, /jl = 0, 6 = 0. Slicing through the graph vertically at a fixed value of y allows us define the fitness 


orderings in Table 3. 
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Sequence 
■No NC 

■Zg,zG,zg 

zG,Zg,zg 
Zg,zg,zG 
zG,zg,Zg 
zg,Zg,zG 
Bzg,zG,Zg 



Figure 5 Fitness orderings. Fitness orderings showing the plausibility of niche construction on the zombie and gooey traits, across parameter 
space. Plots are for y values of 0.1 (panel a), 0.3 (b), and 0.1 (c) ; if not varying along thex- ory-axis, parameters values are fit = \,$c = 0.5, z = 0.2, 
0 = 0.25, v = 0.14,(5 = 0.05, fi = 0,6 = 0. 



construction may maintain the two traits but is unlikely to 
lead to their acquisition (blue region), as going from zg to 
ZG now requires acquiring both traits simultaneously. 

Discussion 

In this study we have explored the potential for niche 
construction to play a role in the acquisition and main- 
tenance of a significant manipulation of host behaviour 
by members of the Baculoviridae family. Results from 
our modelling suggest that the feedback from manipu- 
lating the host is sufficient to acquire and maintain both 
the behavioural manipulation (zombie) and liquefaction 
(gooey) traits. Niche construction provides a strong expla- 
nation for the acquisition of these traits by virus, because 
the ecological inheritance from manipulating the host 
via the expression of the egt gene changes the selective 
pressures over successive viral generations until lique- 
faction changes from neutral or slightly advantageous to 
highly advantageous (Figure 1A). The resulting fitness 
landscape, which includes epistatic interaction between 
the loci, is an effect of the process of niche construction 
that we model here. This process affects the evolution of 
the virus over subsequent transmission events and sea- 
sonal infections, since the virus receives greater fitness 
benefits as it travels through genotype space (Figure IB) 
by acquiring both traits (and possibly fine-tuning their 
expression). It is also worth noting that the system we 
are modelling is an example of the inceptive perturba- 
tion category of niche construction described in Table 



3.1 of [37], wherein "organisms initiate a change in their 
selective environment by physically modifying their sur- 
roundings". This category also applies to the possible niche 
construction effect of egt in which the deactivation of 
moulting in the larvae changes their growth pattern and 
lifespan and thus modifies the occlusion body yield of the 
cadaver. 

Our model results are sensitive to parameter values that 
reflect ecological forces. For instance, differences in host 
ecology across species such as the choice of host plant and 
density of hosts per plant could affect transmission to con- 
specifics or increase contribution to the reservoir, which 
would change the predictions of our model as to whether 
niche construction would favour ZG genotypes or the way 
in which they might have been acquired. Examples of dif- 
fering host ecology abound; for instance, caterpillars of 
Lymantria dispar feed on leaves near the tops of trees, 
providing boosts to horizontal transmission (increasing 
Pl,zg) and favouring ZG. On the other hand, larvae of crop 
pests like Anticarsia gemmatalis (velvetbean caterpillars) 
that feed on soybean and other legume plants have a dis- 
persed spatial distribution [38] and are active enough to 
jump off the plant if disturbed. This may not favour the 
G allele, which may account for the lack of liquefaction in 
this NPV [39], and whether the virus manipulates its host 
is unknown. 

We argue that niche construction provides a more com- 
plete view of the benefits of host manipulation to viruses 
by providing a clear link between traits that, together, 



Hamblin and Tanaka BMC Evolutionary Biology 201 3, 13:1 70 
http://www.biomedcentral.eom/1 471 -21 48/1 3/1 70 



Page 9 of 10 



enhance fitness through the action on the host. Alterna- 
tive explanations, such as standard natural selection with 
epistasis between the zombie and gooey traits, are not 
sufficient because they do not account for the fitnesses 
changing over time in this system; the zombie trait cre- 
ates an ecological inheritance that builds up and cannot 
affect the fitness of the gooey trait until the behavioural 
manipulation has run its course. For Baculoviridae, this 
allows us to explain why the virus would manipulate its 
hosts behaviour in such a dramatic fashion; the result- 
ing change in the host creates conditions that modify 
the effect of the liquefaction trait on viral fitness and 
couple the two traits together. The niche construction 
explanation for viral manipulation of hosts can also pro- 
vide testable predictions for empirical work; our model 
predicts that, when niche construction is favoured, dis- 
rupting the link between between the zombie and gooey 
traits by blocking the manipulation inspired by genes like 
egt should render both (nearly) neutral and thus subject 
to silencing or removal. This may also help to explain 
why in vitro cultures of NPVs often lose their egt genes 
quickly and show much greater levels of virulence [18,19], 
and might provide insights into methods for enhancing 
the effectiveness of baculoviruses as biocontrol agents 
for their lepidopteran hosts, which are often agricultural 
pests. We note that the focus of much parasitological lit- 
erature on the evolution of host manipulation by parasites 
has been on extended phenotype theory as a way of view- 
ing the facts' [40], though the structure of our model 
does not depend on the details of the distinction between 
extended phenotype and niche construction [41,42]. 

Empirically, the family Baculoviridae has been exten- 
sively studied on a molecular level, both for its use as a 
vector for protein expression and as a possible biopesti- 
cide. In addition to this, sufficient basic research on the 
properties of its pathology and ecological inheritance is 
available to make studying niche construction in this sys- 
tem possible; all that remains is to put the pieces together 
appropriately. For example, recombinant strains with egt 
and chitinase/cathepsin deletions have been created (sep- 
arately) for LdMNPV, and the effects of the deletion on 
host behaviour / transmission studied (e.g. egt: [13], and 
chitinase: [43]). The next step is to disrupt the niche- 
constructed pathway in an experimental setting. One pos- 
sible experimental design would be to create the four 
genotypes in our model (zg,ZgzG,ZG) and compare their 
effects on transmission in the lab and in the field, using 
the methodology of D'Amico et al. [43] (confinement of 
larval treatment conditions to mesh bags in a naturalistic 
setting) to measure transmission. Replicating the fitness 
orderings from our model would provide support to our 
predictions. It should also be possible to use an ecologi- 
cal manipulation to test this model. If larvae infected with 
ZG (wildtype) viruses are prevented from experiencing 



the behavioral effects of tree-top disease (i.e. climbing to 
the top of the host plant), then over multiple generations 
we expect rapid loss of egt and (or) chitinase or cathepsin 
in treatment strains but not control strains that are free to 
manipulate the host. 

With over 600 viruses described in the literature [44] 
and mixed success implementing baculoviruses as pesti- 
cides to date [45], finding hosts with ecology favouring 
the niche constructive pathway we have identified in this 
study may aid in designing more effective baculovirus bio- 
control agents. For instance, it may be possible to identify 
which hosts will benefit from increased viral transmission 
via this niche constructed pathway and target the engi- 
neering of baculovirus biopesticides accordingly. It may 
also be possible to manipulate the environment of the host 
to work with the virus. For example, if the effect of zom- 
bie gooey phenotypes is enhanced by a clustered spatial 
distribution of host plants, baculoviruses targeted at crop 
pests may be improved by modifying the planting regime 
to match. Further basic research on the effects of niche 
construction in this viral family thus holds the potential 
to benefit our understanding of viral evolution as well as 
enhance practical applications of Baculoviridae, 

Conclusions 

Niche construction provides a new lens through which 
to study viral manipulations of the host. The analysis 
of the mathematical model developed in this study has 
shown that niche construction provides a plausible expla- 
nation for the behavioural manipulation and liquefaction 
of Lepidopteran hosts by baculoviruses. Understanding 
the evolutionary forces behind these manipulations pro- 
vides not only new insights into viral evolution, but can 
form the basis for new empirical work on baculoviruses 
and research on the use of these viruses as biocontrol 
agents. 
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